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We consider a two-component scaling picture for the resistivity of two-dimensional (2D) weakly 
disordered interacting electron systems at low temperature with the aim of describing both the 
vicinity of the bifurcation and the low resistance metallic regime in the same framework. We 
contrast the essential features of one-component and two-component scaling theories. We discuss 
why the conventional lowest order renormalization group equations do not show a bifurcation in 2D, 
and a semi-empirical extension is proposed which does lead to bifurcation. Parameters, including the 
product zu, are determined by least squares fitting to experimental data. An excellent description 
is obtained for the temperature and density dependence of the resistance of silicon close to the 
separatrix. Implications of this two-component scaling picture for a quantum critical point are 
discussed. 

PACS numbers: 71.10.Ca,71.30.+h,73.20.Qt 



I. INTRODUCTION 

The determination of the low temperature properties, 
and ultimately the ground state phase diagram, of two- 
dimensional (2D) systems of strongly interacting elec- 
trons at low carrier density in the presence of low lev- 
els of disorder remains a very challenging problem in 
spite of a great deal of experimental and theoretical ef- 
fort. The longstanding view that all charge carrier states 
become localized at zero temperature (T = 0) in the 
limit of a large 2D system has been called into ques- 
tion by experimental observations of a finite temperature 
"transition" from insulator-like (dp/dT < 0) to metal- 
like (dp/dT > 0) behavior as the electron density is in- 
creased in very high purity semiconductor MOSFETs and 
heterostructuresi*2iii4. This bifurcation of the resistivity 
p(T) into two families of curves occurs at a critical car- 
rier density of the 2D electron (or hole) system. The 
critical carrier density depends on the specific disorder 
characteristics of the given sample. In the vicinity of the 
bifurcation, p(T) has intriguing scaling properties as a 
function of density and temperature. The lower density 
insulating family of curves having dp/dT < can be col- 
lapsed onto a single curve when plotted as a function of 
a scaling variable T (S)/T. The new density dependent 
temperature scale To(<5) has a power law dependence on 
the magnitude of d — (n — n*)/n* where n* is the critical 
density at the bifurcation. Similarly, the metallic family 
of curves at higher density with dp/dT > collapses onto 
another unique curve with a power law dependence on 
T (S)/T with precisely the same critical exponent. There 
is also scaling behavior with respect to density via \8\ and 
applied electric field, with a different characteristic crit- 
ical exponent. A review has been given by Abrahams et 
aU. 

A bifurcation with scaling behavior at finite T is ex- 
pected to be a generic feature associated with a T = 
quantum critical point (QCP)&. Consequently, the 



demonstrated scaling for the resistivity can be taken as 
evidence for a metal- insulator transition (MIT). With 
this assumption, Dobrosavljcvic et al have given a phe- 
nomenological description of the observed scaling prop- 
erties and power laws in terms of a one-parameter scaling 
theorji In this picture, the bifurcation point lies on a 
separatrix which terminates at T = at a QCP which 
is a fixed point of a set of scaling equations. The fixed 
point is repulsive with respect to the temperature vari- 
able and separates the insulating and metallic regions of 
the low T phase diagram. There is a unique correspon- 
dence between the QCP and the bifurcation point and a 
two-phase ground state is implied. 

However, this correspondence provides only indirect 
evidence for a MIT. Of course the power law fits near 
the bifurcation point are based on linearization of a scal- 
ing flow equation and the linearization must eventually 
fail at sufficiently low temperature, for which Tq(5)/T is 
no longer small, but this is not evidence against a QCP. A 
more crucial point is that a physical effect which is neg- 
ligible for temperatures T ~ 1 K near the bifurcation, 
may become dominant in the low T limit. Scaling prop- 
erties characteristic of an imminent QCP would then still 
be observed near the bifurcation but the suggested QCP 
itself could be totally removed in the zero temperature 
limit. 

A possible example of this situation is the proposal by 
Simmons et a& that in p-type GaAs the p(T) curves of 
the metallic family will eventually turn upward to exhibit 
insulator-like behavior with dp/dT < if the electron 
temperature T could be made sufficiently low. A similar 
conclusion has recently been reached for Si MOSFETs; 9 
However the turn-up is actually only directly observed 
for metallic curves just above the critical density, and, 
as we discuss below, the magnetoresistance data cited as 
evidence of a universal weak localization in 2D in Refs. 
and may alternatively reflect different magnetic field 
dependencies of the competing localizing and dclocalizing 
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triplet spin state interactions. Irrespective of the mecha- 
nism responsible for the upturn, properties of the system 
at temperatures well below the bifurcation temperature 
are not correctly described by a one-component scaling 
theory. We will outline a semi-empirical two-component 
scaling theory that describes within a single framework 
the temperature dependent resistivity both in the vicin- 
ity of the bifurcation and in the low resistance metallic 
regime. 

Following early work by Altshulcr et al*& further 
progress toward an understanding of disordered 2D elec- 
tron systems in the strong coupling regime was made 
by FinkelsteiniiiiS and by Castellani et who showed 
that electron-electron interactions can lead to dp/dT > 
(metallic behavior) in the disordered 2D system at finite 
T. This is in contrast to the localizing insulating behavior 
found in the absence of electron-electron interactions^. 
The physical properties of the system at low T are deter- 
mined by nonlinear interactions of low energy diffusive 
modes. A set of rcnormalization group (RG) equations 
was generated by means of a formal perturbation expan- 
sion in powers of a dimensionless temperature dependent 
resistivity TZ which is related to the physical resistance 
per square by TZ = (e 2 /nh)B^. These scaling results in- 
clude electron-electron interactions to all orders and are 
valid in the diffusive regime fegT < K/t provided TZ is 
small (in principle, 1Z << 1). Additional progress was 
made by Zala et who considered interaction correc- 
tions to transport properties at intermediate tempera- 
tures in the ballistic regime ksT > K/t to all orders in 
the electron-electron interaction using a Fermi-liquid ap- 
proach. These results are also restricted to the perturba- 
tive regime, TZ « 1. Comparisons with experiment for 
the temperature dependent resistivity of Si MOSFETs 
have been made for both the scaling theory^ and the 
intermediate temperature theoryAiAS. The intermediate 
temperature theory gives consistent semi-quantitative 
agreement with experiment in the range of densities for 
which 1Z is small. 

We emphasize that both these theories are perturba- 
tive in 1Z and arc limited to the regime of small 1Z. Hence 
they must fail to describe the regime near bifurcation 
where TZ ~ 1 and they contain no bifurcation point. In 
addition, in the RG procedure the equations develop sin- 
gularities at a nonzero temperature so that the approach 
to the ground state cannot be described. 

In contrast to the 2D case there is a qualitative under- 
standing of the metal-insulator transition in 3D electron 
systems. A consistent qualitative picture of the transition 
has been obtained by RG methods in a space of dimension 
2 + e with e > taken as an expansion parameter. Physi- 
cal quantities such as critical exponents are expressed as 
power scries in e. While only the leading terms in e have 
been obtained it is presumed that resummation meth- 
ods would give good results in 3D if sufficient correction 
terms could be calculated. In the case of e > it is clear 
that the interplay of strong electron-electron interactions 
with disorder is an essential aspect of the bifurcation. 



It is remarkable that in spite of great effort a cor- 
responding theory does not exist for 2D systems, that 
is for e = 0. This lack casts doubt on the physical 
significance and relevance of the "universal" aspects of 
electron-electron interactions in disordered systems. In 
addition the lack of even an approximate scaling theory 
has prevented quantitative discussion of the low temper- 
ature limit and the approach to the ground state. We 
expect that the interplay of strong electron-electron in- 
teractions with disorder must be an essential aspect of 
the bifurcation in 2D irrespective of whether or not a 
true metal-insulator transition occurs. Of course in some 
systems material dependent effects also contribute. 

In this paper we discuss why the lowest order RG equa- 
tions in 2D do not show a bifurcation and we give a 
semi-empirical extension which does describe a bifurca- 
tion region. We focus on the immediate vicinity of the 
bifurcation and on the metallic regime. The strong in- 
sulator limit will be discussed elsewhere. The proposed 
scaling equations have a physical low temperature limit 
and have no singularities at finite length or temperature 
scales. This scaling picture provides a theoretical frame- 
work for the interpretation of experimental results, more 
specifically the "universal" contributions due to electron- 
electron interactions and disorder. 



II. RENORMALIZATION GROUP EQUATIONS 

The RG equations established at one-loop level for a 
2D system of electrons in the presence of disorder are 
based on four dimensionless scaling parameters TZ, 72, 
Z and 7 c( iiiiSii^ where T2 = Z72 is the electron-hole 
scattering amplitude for the triplet spin state, r c = Z^ c 
is the singlet state particle-particle scattering amplitude 
and Z is the dynamical energy rescaling function. To- 
gether with TZ, these quantities are all functions of the 
variable y = In A -1 which describes rescaling of the mo- 
menta after integrating over the momentum shell speci- 
fied by Afcg < k 2 < k 2 m*L We consider only the case when 
the disorder is due to purely potential scattering. When 
a connection between length and temperature scaling is 
needed, we relate a thermal length £ t h to temperature 
by W4z = (Tei/T) 1 /*, where t el and T el are the length 
and temperature scales for elastic scattering and z is the 
dynamical critical exponent^. 

Particle-particle scattering is not considered to play an 
essential role relative to the electron-hole scattering rep- 
resented by 72 so we omit consideration of "f c . The energy 
scaling function Z(y) is not followed explicitly but its ef- 
fect is taken into account at the bifurcation fixed point 
by y = ln^Tej/T) 1 / 2 ]. These approximations permit dis- 
cussion of a two-component scaling theory based on the 
scaling parameters TZ and 72. With these simplifications 
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the RG equations are 

dTZ/dy 
d-y 2 /dy 

where 
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n v is the number of valleys, so for Si n v = 2. Equations 

(1) and (2) arc to be integrated upward in y, starting 
from initial bare values for the dependent variables at 
y = y(0) for which T e i/T is of order unity. Increasing 
y corresponds to integrating out shorter wavelength and 
higher energy excitations. This amounts to increasing the 
length scale and decreasing the temperature scale. Our 
Eqs. (1) - (3) are essentially the same as Eqs. (1) and 

(2) of Ref. [Hi The notational difference is that we use a 
running variable y = — ln(A) as in Ref. H^. while Ref. ^| 
uses £ = — ln(fcsTr/?i). Prior to renormalization due to 
interacting diffusive modes, we can set z = 1. After the 
RG flow has proceeded to the bifurcation fixed point, we 
set y = IniiTei/T) 1 ^} where z is the dynamical critical 
exponent of the fixed point. 

The first term (unity) on the right hand side of Eq. 
Q contains no electron-electron interactions and arises 
from weak localization. On scaling this term enhances 
the resistivity. The other term, in the square brackets, 
is due to electron-electron interactions in the singlet and 
triplet particle- hole spin states and can be negative when 
72 is large enough. When it is negative this term has the 
opposite trend of reducing the resistivity upon scaling. 
The bare value of 72 is given by ^ are = — Fq/(1 + Fq) 
where Fq < is the spin-antisymmetric Landau param- 
eter. From Eq. (2), 72 increases with rescaling. With 
n v = 2, the expression for 0(72) changes sign when 72 
reaches 0.46. 

A change in sign of a (72), corresponding to a net delo- 
calizing effect in zero external magnetic field, has no par- 
ticular a priori implication for magnetoresistancc. The 
magnetoresistancc data of Refs. |8||9| can be consistent 
with our scaling conclusions if, for example, the localiz- 
ing contribution in Eq. @ has a stronger magnetic field 
dependence (at small fields) than the delocalizing contri- 
bution from the electron-electron interactions. Thus the 
appearance of a weak localization precursor signature in 
magnetoresistancc data does not necessarily imply there 
will be an eventual turn-up in the zero field resistivity at 
very low temperatures. 

As we integrate Eq. (2), 72 increases from its initial 
bare value but diverges at a finite value of y = y ma x 
provided 1Z remains finite at y = y max Mii^^ It is easy 
to confirm that a solution of Eqs. (1) and (2) with finite 
1Z = TZ m in and 72 arbitrarily large is consistent. In this 
limit Eq. (2) can be written 



It follows that 72 diverges at a finite y m ax- To verify the 
consistency of a finite lZ m in when 72 diverges, we can 
divide Eq. (1) by Eq. (2) and rearrange obtaining, 



dK/K= -d72(61n 72 )/72 2 



(5) 



dj 2 /dy = (JZ ri 
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The integral of the left hand side is finite so it follows that 
1Z indeed reaches a finite lower limit lZ m i n as 72 — > 00. 
The rescaling cannot be continued beyond the singularity 
at y = y m ax which means that the zero temperature limit 
cannot be reached. 

In principle, this singularity in the triplet state scatter- 
ing amplitude might signal the onset of a magnetic insta- 
bility in the system. On the other hand, the divergence 
might simply be an artefact of a low order perturbation 
expansion. The ultimate fate of this singularity at very 
low T is not completely clear. Castcllani et showed 
that a similar divergence in the energy rescaling function 
Z(y) can cause the singular point to shift to extremely 
low T. The paramagnetic metallic regime would then 
extend over a wide temperature. More recently, Kirk- 
patrick and Belita 2 ^ and Chamon and Mucciolo^i have 
indeed found a solution of the RG equations correspond- 
ing to a disordered ferromagnet. The question of how 
low the temperatures would be where such transitions 
might occur is open. In the absence of theoretical guid- 
ance on this point we turn to experiment for information 
on the triplet spin state scattering amplitude 72 in the 
temperature range of the bifurcation. 

References concluded from an analysis of magne- 
toresistance data that there would be a ferromagnetic 
instability very near the densit y of the bifurcation. How- 
ever, measurements in Refs. |23j|24j found an enhancement 
of the effective g- factor g* , but no singularity for densities 
down to r s = 8.4, a range which includes the bifurcation. 
This is interpreted in terms of the spin-antisymmetric 
Landau parameter as g* — 2/(1 + Fq ). Fq would be — 1 
at the onset of a ferromagnetic instability. For r s = 8 
Ref. I24I gives Fq ~ —0.5, corresponding to a spin suscep- 
tibility enhancement factor of 2. At temperatures well 
above the bifurcation the bare values of the electron- 
electron interaction amplitudes like ^ are have negligi- 
ble diffusion corrections. 7^ are is then related to F a by 
7 &are = -F§ / {1 + Fq 2 ). A value Fq = -0.5 corresponds 
to 7& are = 1. On the basis of the direct measurement 
of g* we will assume that the triplet spin state scatter- 
ing amplitude 72 is finite and well behaved throughout 
the density and temperature range of the experiments we 
consider. With this assumption, an explicit RG equation 
is not needed for 72 . Nonmagnetic Fermi liquid behavior 
of the system is a sufficient condition for a smooth 72 but 
may not be necessary. 



III. BIFURCATION IN 2D 

To discuss why Eqs. (1) and (2) fail to describe a bi- 
furcation in 2D electron systems, we first recall the corre- 
sponding RG results in d = 2 + e dimensions with e small 
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and positive. Making the same physical assumptions re- 
garding the interacting diffusive modes and in the same 
one- loop approximation. Eq. (1) becomes 

dTZ/dy = -{e/2)TL + a{ l2 )TZ 2 . (6) 

The first term on the right hand side of Eq. JHJ) is the 
consequence of the L 2 ~ d form factor when converting re- 
sistivity to resistance in a space of d dimensions, and the 
coefficient 0(72) in the second term is the same as in Eq. 

We identify the bifurcation point as the point at which 
dTZ/dy = 0. The zero of Eq. © occurs at the criti- 
cal value TZ* = e/(2a), provided a > 0, where a is the 
value of a (72) when the bifurcation occurs. In order for 
the bifurcation to be a precursor for a quantum critical 
point, the temperature must be a relevant variable. This 
requires that t _1 > in the linearized flow equation 

d{K-K*)/dy = t- 1 (K-K*) . (7) 

This procedure for identifying a critical point by a 
linearized flow equation is standard. From Eq. JHJ, 
r e/2 > 0. The condition that t _1 > implies that 
the resistance has positive (negative) slope with respect 
to T in the metallic (insulating) regime. 

It is clear that a non-trivial root of dTZ/dy = re- 
quires at least two terms and that for the 2D case Eq. 
(1) as presented (that is, Eq. © with e = 0) will not be 
sufficient. Attempts have been made to discuss possible 
metal-insulator transitions in 2D using Eq. (1) by tuning 
the parameter 0(72) to zero. However, a critical point 
is a robust property of the entire dTZ/dy and cannot be 
described by the properties of only a single term a(j 2 )TZ 2 
in a series. Such a procedure is not stable to the addition 
of higher order terms. The simplest modification of Eq. 
(1) that can show a bifurcation and is consistent with 
two-component scaling is 

dTZ/dy = a(j 2 )TZ 2 + p{ l2 )TZ 3 + ... (8) 

The function ^(72) is not known explicitly but its sign 
can be determined by the conditions that TZ* and r _1 
are both positive. Linearizing Eq. (jSJ) about the zero of 
its right hand side gives the linearized flow equation Q 
with TZ* — — a 1/5 and exponent r = (3/ a 2 . Here a and 
(3 are the values of 01(72) and ^(72) when the bifurcation 
occurs. TZ* and t _1 are both positive provided a < and 
(3 > 0. There is no bifurcation in 2D if 0(72) is always 
positive. Integrating Eq. Q starting from an initial TZo 
at y = yo gives two families of curves. A metallic regime 
and a bifurcation can thus be described. 

We conclude that the bifurcations in the 2D system 
and the 2 + e system are controlled by different fixed 
points. The fixed point TZ* = e/a for 2 + e requires a > 
which is in the range of weak electron-electron interac- 
tions. This fixed point becomes trivial (TZ* = 0) in the 
e — > limit and plays no role in 2D. This allows a new 
fixed point TZ* = —a/f3 to become physical in 2D at 
a scale where electron-electron interactions have become 



strong enough to change the sign of a. In each case there 
is only one fixed point and the physical picture of inter- 
acting diffusive modes is correct for the determination of 
the universal contributions. 

The strength of the electron-electron interactions in 2D 
is crucial for generating a bifurcation. At high densities 
where the electron-electron interactions are weak, j^ are is 
small and a(j 2 are ) is positive. As the density is lowered 
the initial ^ are increases. Using the data of Ref. f° r 
g* in Si, the sign change in a(j 2 are ) occurs for a density 
corresponding to r s ~ 3. This provides an upper limit to 
the density at which a bifurcation can occur in Si. For 
r s < 3 the electron-electron interactions are too weak. 

In order to describe quantitatively the resistivity near 
the bifurcation where R ~ 1, as well as in the insulating 
regime where R >> 1, the sum of the series implied in 
Eq. 1JSJ must be adequately represented. If the series is 
truncated at an arbitrary finite order a spurious diver- 
gence at a finite y m ax can occur in the insulating region. 
While in the insulating limit TZ and dTZ/dy are expected 
to diverge as y — > 00 (that is, at T — > 0), the divergence 
in dTZ/dy must be sufficiently weak that a spurious diver- 
gence in TZ at a finite y m ax does not occur. A linear power 
law in TZ is the strongest growth of dTZ/dy at large TZ that 
permits this (with possible logarithmic corrections). For 
simplicity, we maintain the fixed point structure of the 
low order terms and introduce a denominator into Eq. (jHJ) 
to represent the net effect of higher order terms including 
the linear growth at large TZ, 

d/TZ a( l2 )TZ 2 + P( l2 )TZ 3 

dy l + K ( l2 )TZ 2 ■ [) 

Equation ||§J) satisfies the minimal conditions of having 
a bifurcation with two classes of well defined solutions 
(metallic and insulating) depending on the choice of ini- 
tial conditions at y = y(0). It may be regarded as a 
semi-empirical Padc approximation to the full series. Of 
course, these minimal conditions do not uniquely deter- 
mine the functional form (see also Section IV. B). 

We focus in this paper on the metallic regime and the 
close vicinity of the separatrix. The known perturba- 
tive results are contained explicitly in Eq. © so the 
solution for y — > 00 and TZ — > is exact. The strongly 
insulating limit TZ — » 00 contains additional logarithmic 
corrections so Eq. @ is incomplete in this limit, and a 
detailed discussion of the deep insulating regime will be 
given elsewhere. 

Close to the separatrix the functions 0(72), (3(j2), and 
^(72) are taken to be slowly varying and so are replaced 
by their constant values a, (3 and k for 72 near the start of 
the bifurcation. Information on these parameters is given 
by fitting to experimental data in the following section. 
Equation @ can then be rewritten in the form 

1 dTZ _ I A + A 2 

TZ~dy~ ~ 7l + 0(2A + A 2 ) ' (1 ' 
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where 

TZ* = -a/0 

tz/tv = 1 + A 

t- 1 = (3TZ* 2 /{1 + kTZ* 2 ) 

<j) = kK* 2 /(1 + kK* 2 ) . (11) 

Linearizing Eq. (|10f> in A we recover Eq. Q, with the 
solution 

\n{\TZ - K*\/\K(0) - TV\} = r" 1 ^ - y(0)) . (12) 

Since r describes the rescaling of an inverse length 
squared, it is related to the critical exponent v of the 
correlation length by r = 2v. The temperature is intro- 
duced by the thermal length tth ~ l/T 1 / 2 giving 

\K-K*\ = \K(0) -K*\ (T /T) 1/Zl/ , (13) 

with the same exponent for both the metallic and the 
insulating branches. The prcfactor defines a tempera- 
ture scale T (S) ~ \5\ zv T a . Both of these features agree 
with the observed scaling and with the phenomenological 
scaling of Ref . 0. 

IV. RESULTS 
A. Vicinity of separatrix 

We now compare our results based on Eq. (|1U|) with 
experimental values of the resistivity p(T) close to the 
separatrix obtained for a Si-MOSFET taken from Fig. 
1 (b) of Ref. [25| (see also Ref. @). Equation O can be 
integrated analytically. The parameters TZ* , z, z/, a, (3, 
and k in TZ(y) are chosen to give a best fit to the exper- 
imental data. Not all these parameters are independent 
of each other. The experimental data fix TZ* = 2.8. Since 
TZ*=—a/f3, we can then consider a to be fixed with /3 the 
free variable. The z and v = t/2 enter together as a prod- 
uct via T~ x {y- y(0) = ln(T /T) 1 /* i/ . The temperature T 
just prior to the bifurcation was taken to correspond to 
the temperature scale of elastic scattering T e i = 1.75 K. 
From the definition of r, we have k = t(3 — 1/TZ* 2 : which 
relates k to /3. There are therefore three independent 
variables z, v = t/2 and (3. 

From the combination of electric field scaling and tem- 
perature scaling of the resistivity the dynamical critical 
exponent z is believed to be in the range 0.8 to 1.2£ We 
have therefore made least squares fits to the experimental 
data with b = zv and j3 as free parameters for the fixed 
values of z in the range 0.8 < z < 1.2. The optimum 
values of the fitting parameters are determined by mini- 
mizing the root mean square relative deviation between 
theory and experiment 

fl = Iyv!2 *j ' ( i4) 

3=1 3 



where N is the total number of points included in the fit. 

We have also examined the sensitivity of the fitted pa- 
rameters to the temperature and density range of the fits. 
This is essential because the values of the least squares 
fitted parameters can vary with the range of temperature 
and density considered. Eq. (jlO|l has been established 
on the assumption of constant parameters for sufficiently 
small Tq(S)/T. Fits to data will be valid only if the de- 
rived parameters are stable with respect to reducing the 
maximum allowed Tq(S)/T, that is by restricting the tem- 
perature and density range of the fits. 

We first carried out fits including all of the data points 
in Fig. 1(b) of Ref. l25l To test sensitivity to the range 
of fit we successively restricted the allowed data points 
by the conditions T (5)/T < 1.0, 0.5, and 0.25. For 
Tq(S)/T < 0.5 and fixed z the D and the values of the 
fitted parameters become stable. The coefficients (3 and 
k vary slowly with z. However the essential parameters 
from the point of view of fitting to a universal scaling 
form are constant throughout the range of z. These val- 
ues are D = 0.036 ± 0.0005, b = zv = 1.09 ± 0.005 and 
<f) = 0.80±0.005, where the uncertainties reflect the small 
variations due to z. A constant b and 4> is consistent, 
since (2(f) — l)/b is the universal coefficient of the first 
non-linear correction A 2 in Eq. (HOC . 




0.5 1.0 1.5 

T(K) 



FIG. 1: Solid lines: measured resistivity p(T) in units of h/e 2 
from Ref. [25] in a Si-MOSFET close to separatrix for electron 
densities (from the top) n = 8.6, 8.8, 9.0, 9.3, 9.5, 9.9, 11.0 x 
10 10 cm~ 2 . Dashed lines: our calculations using Eq. (10). 
The values of the parameters are given in the text. Dotted 
lines are in the small T region Tq(8)/T > 0.5 that is excluded 
from the fit. Small vertical arrows: see text. 

Figure 1 compares our results (dashed lines) from 
Eq. (|10f) for the dimensionless resistivity in the form 
pe 2 /h = TZir with the experimental measurements (solid 
lines). The agreement is excellent. The dotted lines 
indicate the data points excluded from the fit by the 
condition Tq(S)/T < 0.5, and the small arrows show 
the edge of the corresponding excluded data range for 
T (5)/T < 0.25. 
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B. Exponential form in metallic region 

Previous fits to experimental data have shown that 
the temperature dependence of the resistivity is approxi- 
mately exponentialSSiSl. This result was accounted for by 
Ref. 13 by arguing that the beta-function of the conduc- 
tivity was a logarithmic function of conductivity, even in 
the metallic range near the separatrix. This is equivalent 
to assuming that the sum of the full series in Eq. JHJ, 
again with the coefficients evaluated at the bifurcation, 
is approximately logarithmic. Then Eq. (|10|) is replaced 
by 



giving 



dlnTl/n* 
dy 



Tl{T) ~ TV exp 



r" 1 In ft/ft* , 





\T (sy 






T 





(15) 



(16) 



with C a positive constant of order unity. 

For the purposes of least squares fitting to data for 
small |A|, it is important to note that Eqs. IjlUf) and 
(|15f) are not inconsistent. Standard procedures for Pade 
approximations give the identity 



ln(l + A) 



1 + A- 



C(A 5 ) 



(17) 



The exact sum of the perturbation series near the bifur- 
cation has the same form 



d\nTl/TV 
dy 



1 



aA 2 



t 1 + bA + cA 2 



0(A 5 ) 



(18) 



with the coefficients being of order unity. Equation 11U|) 
and Eq. (fTST) with (|T7|) can thus be regarded as different 
approximations to Eq. (|18|) and are therefore equivalent 
from the point of view of least squares fitting to data for 
small |A|. It is interesting to note that the low order 
expansion of Eq. JT7J), ln(l + A) = A - ±A 2 + C(A 3 ), 
gives a value (f> = | which is comparable to our fitted 
value of 6 = 0.80. 



C. Deep metallic limit 

Due to its y dependence 72 continues to rescale as 
the low temperature limit of the metallic regime is ap- 
proached. If we assume that 72 approaches a finite limit- 
ing plateau value 7° as T — > 0, then 0^(72) and ^(72) also 
approach finite limiting values 01(72) and flij^)- Simi- 
larly, again due to the rescaling, the zero temperature 
limit z°v° of the exponent zv is expected to be different 
from the bare value and also from the value determined 
at the bifurcation. In this case for T sufficiently small we 
can obtain a consistent solution of Eq. JHJl with TZ << 1 



so that the perturbation expansion is well represented by 
the leading term, 

dTZ/dy = a(j%)1l 2 + 0{TV) , (19) 

with the solution 

K/K = [l-a(72°)(2/-2/o)] _1 • (20) 



Using (y-yo)-- 
conductivity g 



h^Tb/T) 1 / 2 and with the dimensionless 
= TV 1 , we find 



dg 

d\nT 



Mil) 



(21) 



The dominant temperature dependence is then logarith- 
mic as observed experimentally. References |53 and 
express this logarithmic contribution to the conductivity 
G as AG = (e 2 /h)C{n)\nT. Using Eq. JSTJ we identify 
C{n) = a(72)/(7rz°). As the density increases 7° and 
a (7°) decrease in magnitude and 0(72) may even change 
sign. This dependence of C{n) on density is in agreement 
with that observed experimentally. 



V. CONCLUSIONS 

The low order perturbative RG equations of Refs. 0, 
do not describe a bifurcation in 2D, but a proper 
description of a fixed point and a bifurcation in 2D can 
be obtained when higher order terms in the perturbation 
expansion for 7Z are retained (see Eqs. ©). The 
resulting 2D fixed point exists only for 0^(72) negative 
and so is unrelated to the fixed point in dimensions 2 + e 
for which a(72) is positive. 

Our results provide a coherent semi-empirical two- 
component scaling description of the density and temper- 
ature regime near the observed bifurcation and through- 
out the metallic regime. Near the separatrix a least 
squares fit to experimental data using Eq. IjlOjl gives an 
excellent description of the observed density and temper- 
ature dependence of the resistivity. A full discussion for 
the insulating range will be given elsewhere. The scal- 
ing results apply only to the "universal" contributions to 
the resistivity which are a generic consequences of the 
interplay between electron-electron interactions and dis- 
order. There are also "non-universal" contributions to 
the resistivity which will be material-dependent. 

The present picture is based on a K which shows a 
bifurcation and on a 72 which scales smoothly with tem- 
perature. The bifurcation "point" TV = — a/ (3 therefore 
varies smoothly with temperature and the separatrix, at 
which dTZ/dT = 0, is "tilted" upwards. A tilted sep- 
aratrix can lead in the metallic regime to a turn-up of 
the resistance at low temperatures. This cannot occur in 
a one-component scaling theory where the separatrix is 
flat and the bifurcation point is uniquely determined as 
a function of density. 

An exact solution of Eq. i|19fl in the very low tempera- 
ture limit shows a logarithmic dependence of the conduc- 
tance on temperature with a coefficient which is negative 
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and decreases in magnitude as the density increases. The 
low temperature behavior is in agreement with experi- 
mental results^LSS. The physical origin of this logarith- 
mic electron-electron contribution to the conductivity is 
the same as that observed at high temperature. However 
the numerical value of the low T coefficient can differ 
from the bare value at higher temperature. 
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